Giant Kohn anomaly and the phase transition in charge density wave ZrTes 
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A strong Kohn anomaly in ZrTes is identified in the mostly transverse acoustic phonon branch 
along the modulation vector qp with polarization along the a* direction. This soft mode freezes to 
zero frequency at the transition temperature Tp and the temperature dependence of the frequency 
is strongly affected by fluctuation effects. Diffuse x-ray scattering of the incommensurate super- 
structure shows a power law scaling of the intensity and the correlation length that is compatible 
with an order parameter of dimension n = 2. 



Dynamical modulations of the spin or charge density 
are a ubiquitous phenomenon in crystalline metals and 
form an important ingredient in solid state theory [Ij. 
In many materials they remain dynamical at all tem- 
peratures and contribute to the rich phenomenology of 
electronic, magnetic and thermal transport and order- 
ing. Recently, their role has been re-discussed in high Tc 
super-conducting and colossal magneto-resistance oxides, 
in connection with the formation of striped charge order 
[2]. Similar to the stripes, a static charge density wave 
(CDW) removes spectral weight from the Fermi level, 
while the superconductivity can be supported by dynami- 
cal density waves. In general, the modulation wave vector 
q can be incommensurate with the underlying lattice and 
is derived from a nesting geometry of the Fermi surface 
that provides strong screening at q = 2kp. Such a nested 
Fermi surface is found in particular in low-dimensional 
crystal structures with chain-like atomic arrangements 

Due to the low-dimensional nature, the phase tran- 
sition, the so-called Peierls transition, to the statically 
modulated ground state of the one-dimensional chain is 
strongly affected by fluctuations [3]. Inelastic neutron 
scattering studies revealed a giant Kohn anomaly in KCP 
[3] and blue bronze [5j, whereas (TaSe4)2l [HIT] displayed 
a temperature dependent soft mode, which did, how- 
ever, not attain zero frequency at the Peierls transition 
temperature Tp. An inelastic x-ray scattering study of 
the prototypical ID compound NbSea failed to show any 
clear Kohn anomaly but only a broadening of an acoustic 
phonon mode at the first modulation vector qi of NbSea 

The modulated ground state itself is readily observed 
by the appearance of new diffraction peaks, with a po- 
sition in the reciprocal lattice that represents the mod- 
ulation vector q. Similar to other structural phase tran- 
sitions, the temperature dependence of these superstruc- 
ture reflections is used to study the ordering of the CDW 
modulation. Again, only a small number of ID inorganic 
chain systems show this transition and very few studies 
exist. These find that the Peierls transition is governed 
by a critical exponent of the order parameter, which is 



compatible with an order parameter of dimension n ~ 2 
(blue bronze 9 ) or close to being a (blurred) first order 
transition ((TaSe4)2l [10]). A clear assignment to a par- 
ticular universality class has, however, not been found 
and the question remains open whether the transition 
can be classified in this way. 

While the CDW ground state and its many fascinat- 
ing properties are well described by theory [I] , the phase 
transition is accessible only by experiment. In this letter 
we present an experimental study of the fiuctuations and 
the soft mode phonon in ZrTca to elucidate the proper- 
ties of the charge density wave in the important regime 
close to and above the phase transition, where the density 
waves are of dynamical nature. 

ZrTe3 crystallizes in a monoclinic structure with space- 
group P2i/rn (lattice parameters a — 5.89 A, b — 3.93 A, 
c = 10.09 A, a = 7 = 90°, /3 = 97.8°). (ZrTe3)oo 
chains similar to NbSe3 run along b. Only one type of 
chain is present in the unit cell. The Peierls modula- 
tion qp = (0.07, 0, 0.333) (reciprocal lattice units) [TT] 
is transverse to these prismatic chains with a small a* 
component and a tripling of the unit cell along the lay- 
ering c* direction. The modulation is thus very different 
from NbSe3, and mostly in the almost equidistant Te-Te 
chains along a, where the prismatic chains are laterally 
bonded. Resisitvity measurements show an anomaly due 
to the transition at Tp = 63 K and the system remains 
metallic below Tp [12]. Electronic structure calculations 
[131 E] and a detailed photo-emission study [TS] explain 
this behaviour by nesting in a small electron pocket of 
the highly directional Te 5px band running along these 
Te-Te chains, while other sheets of the Fermi surface re- 
main unaffected by the transition. The opening of a gap 
in the electronic dispersion of this band was found to fol- 
low a BCS model of the Peierls transition with a mean 
field transition temperature Tmf about four times higher 
than Tp |15| . Last, but not least, the phase transition 
gives rise to a very small specific heat anomaly at Tp 

m- 

The experiment was performed at the inelastic x-ray 
scattering (IXS) beamhne ID28 at the European Syn- 
chrotron Radiation Facility. Phonon dispersions were 
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FIG. 1: Acoustic phonon along qp at T = 100 K. (a) Rep- 
resentative IXS spectra at momenta close to (401) at qp 
(a* = 3.93) and across qp on a logarithmic scale with ar- 
bitrary offset. Error bars are smaller than the symbols ex- 
cept for the outer data points. Solid lines below that data 
are (shifted) resolution functions that were summed to fit the 
data, (b) Phonon dispersion extracted from fits to the IXS 
data. The dashed line is a sinusoidal model dispersion Uni {q) 
that serves as a guide to the eye. The inset shows a portion of 
reciprocal space with the scan direction through qp indicated. 

measured with an energy resolution 5E =1.7 meV. The 
momentum resolution was SQ = 0.05 A in the hor- 
izontal a* — b* plane and 0.15 A in the vertical c* 
direction. Diffuse scattering and diffraction scans were 
recorded with an energy-integrating detector at the same 
wavelength of A = 0.57 A and similar momentum resolu- 
tion. The samples were grown by chemical vapour trans- 
port and formed platelets of 1 x 2 x 0.1 mm"^ along a, b and 
c respectively, that were glued to the copper cold finger 
of a closed cycle cryostat using silver paste. The abso- 
lute temperature scale was calibrated relative to a diode 
thermometer mounted on the cold finger by setting the 
inflection point of the intensity of critical scattering at 
Tp to the known Tp = 63 K (see Fig. [I]). 

The phonon dispersion along the CDW modulation 
vector was determined as shown in Fig. [T] IXS spec- 
tra were recorded for each momentum transfer Q close 
to the strong Bragg spot G = (401), along a straight 
line from (401) through qp and beyond, as shown in the 
inset of Fig. [ij In the figure, the a* component of Q 
is indicated in relative lattice units. At a temperature 




3.88 3.90 3.92 3.94 3.96 3.98 4.00 
along c/p (a* component) 



FIG. 2: Phonon dispersion from (401) in the direction of qp 
at various temperatures T > Tp above the transition. The 
inset shows the phonon energy huq at (fp as a function of 
reduced temperature (T — Tp)/Tp. The solid and dashed 
lines represent a 1/8 and a 1/2 power law, respectively. 

T — 100 K = 1.6-Tp, the spectra are dominated by 
a central peak at zero energy loss for all momenta, in 
particular close to qp [a* ~ 3.93) and on approaching 
the F-point (lowest curve). The spectra are fit by a sum 
of three or five terms for one or two phonons, respec- 
tively. The first term is the central peak, the following 
terms are phonons in Stokes- and anti-Stokes positions 
with their intensity ratio given by the detailed balance 
at the corresponding temperature. No convolution was 
performed, and the model function was a sum of shifted 
resolution functions according to these terms, thereby 
assuming that broadening effects will be small with re- 
spect to the experimental resolution, as they cannot be 
resolved in the strongly overlapping peaks. This method 
yielded smooth and self-consistent phonon dispersions. 

The thus extracted dispersion is displayed in Fig.[lj3). 
It shows a dip &t qp deviating from a sinusoidal model 
dispersion. This Kohn anomaly (KA) is a direct manifes- 
tation of the Fermi surface nesting vector 2kp. No such 
anomaly was observed for other acoustic phonons; in par- 
ticular a scan along the high symmetry direction a* (lon- 
gitudinal acoustic phonon) did not show any anomaly, 
although the Fermi surface could also afford a certain 
nesting along this direction. All other modes that were 
studied did not show any anomalies, and also a Raman 
scattering study of the phonon modes at the F-point did 
not show any anomalous temperature dependence [T7]. 
The line along qp corresponds to an almost transverse 
dispersion along the large c* component. The phonon 
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FIG. 3: Diffuse scattering along qp (a) Symbols represent the intensities extracted from data fitting of the IXS spectra: central 
peak (triangles), phonon contribution (filled circles). The solid line is the total diffuse scattering intensity as measured with an 
energy integrating detector, (b) Diffuse scattering intensities at three selected temperatures as indicated, (c) Scans of diffuse 
scattering across qp along b* at various temperatures above Tp. 



polarization is along a*. Thus, the displacement pattern 
of the soft-mode is of mostly transverse character with 
respect to the Te-Te chains whose Te 5px orbit als form 
the quasi 1-dimensional Fermi surface. 

The temperature dependence of the Kohn anomaly is 
traced in Fig. [2j As mentioned above, already at room 
temperature (4.6 • Tp) a slight dip is visible, barely out- 
side the experimental error bars. As the temperature is 
lowered, the dip gets more pronounced, until the phonon 
frequency drops rapidly to zero at Tp. Mean field the- 
ory predicts this freezing of the phonon as a logarith- 
mic divergence that would correspond to a square-root 
power law close to the phase transition [T]. Inspection 
of our data by plotting on a double-logarithmic scale did 
not reveal a well-defined power-law behaviour. Over a 
broad range in temperature, the closest correspondence 
to a power law approximately (T — Tp)^/^. On the lin- 
ear scale in Fig. [2j this 1/8 power law and the theoretical 
expectation of a 1/2 power law are shown. Thus, we con- 
clude that the soft mode freezes in the fashion of a giant 
Kohn anomaly. Our results suggest that the phonon fre- 
quency drops to zero at the phase transition. Above the 
transition, fluctuation effects lead to a much more rapid 
recovery of the non-interacting phonon frequency than 
the logarithmic dependence predicted by mean field the- 
ory. 

The deviation from mean field behaviour led us to in- 
vestigate the thermodynamic signature of the phase tran- 
sition by the study of the diffuse scattering (T > Tp) 
and superstructure intensity (T < Tp) at qp. Again, the 
experimental difficulty lies in the close proximity to the 
Bragg reflection due to the small component a* of qp. 
A survey of momentum space slightly above Tp showed 
that the diffuse scattering signal is present for Q = G±qp 
when G has a strong component along the a* direction 
of the reciprocal lattice. 

Our combined diffuse and inelastic x-ray scattering 



study allows disentangling the static and dynamical con- 
tributions to the diffuse intensity. On inspection the IXS 
spectra show that the dominant contribution is the cen- 
tral peak (zero energy loss). Figure |3^) compares the 
intensity of the central peak and the phonon (symbols) 
with the energy integrated signal (line) at T = 73 K > 
Tp. The central peak is not only dominant, but its mo- 
mentum dependence corresponds to that of the diffuse 
scattering signal. Thus we can interpret the diffuse scat- 
tering as a quasi-static diffraction feature which is repre- 
sentative of the (fluctuating) order of the CDW. 

Scans of the integrated scattering intensity along qp 
are displayed in fig. [3|d) and show the emergence of a dif- 
fuse scattering feature. Due to the strong background 
of thermal diffuse scattering (TDS) from the acoustic 
phonon (dominating at high temperatures), the peak 
shape is not easily determined. A reliable measurement, 
however, is possible for scans through qp along h* as the 
rocking curve along this direction is sharp and no resolu- 
tion broadening affects the width of the intensity profiles. 
Again, at high temperature a broad peak is partially due 
to TDS, but the emergence of a sharp feature with de- 
creasing width at lower T is due to the longer correlation 
length of the CDW order. Below Tp this peak does not 
sharpen any further, although it remains slightly broad- 
ened with respect to the experimental resolution. Its in- 
tensity grows rapidly on lowering the temperature below 
Tp. 

The temperature evolution of the intensity and the 
width along h* are shown in Fig. |4] The transition ap- 
pears slightly broadened around Tp. Outside of this 
range, the intensity and width show power law scaling 
with reduced temperature. The critical exponents de- 
rived from these are 2/3 — 0.26 ± 0.06 for the intensity, 
that is taken here as a measure of the order parameter 
squared, and v = 0.85±0.2 for the width, that is propor- 
tional to the inverse correlation length. The diffuse scat- 
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tering peaks as well as the diffraction profiles below Tp 
were always well represented by a single Lorentzian func- 
tion. Other critical exponents could not be derived, in 
particular the shape of the diffuse scattering peak along 
a* and c* is not accessible due to its large width on a 
sloping strong background. 

The critical exponent /? = 0.13 ± 0.03 is even smaller 
than the ones observed in other CDW materials (0.32 ± 
0.05 for blue bronze i9j and 0.23±0.02 for (TaSe4)2l llOj. 
Thus, a mean field description of the phase transition is 
not appropriate and fluctuation corrections are needed. 
Such small exponents (3 of the order parameter are best 
compatible with a two-dimensional order parameter. The 
diffuse scattering spot is well defined in momentum space 
and the Kohn Anomaly is sharp in a three dimensions, 
thus above Tp a three-dimensional fluctuation regime is 
observed. Furthermore the exponent v of the correlation 
length is compatible with a three-dimensional XY univer- 
sality class [TS] . On the other hand, the Peierls transition 
could be considered as a discontinuous freezing (1st or- 
der) that is blurred by the presence of defects and thus 
rendered continuous. This would also lead to a power 
law with a small exponent |10j . However, our experi- 
ment did not show any hysteresis in temperature scans 
and only a slight specific heat anomaly was reported in 
ZrTea [T6|. Thus, we conjecture that the phase transition 
is continuous (2nd order) or very close to that. 

In summary, we have identified the soft phonon mode 
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FIG. 4: (a) Temperature dependence of superstructure inten- 
sity at — (left scale) and the width of the diffuse scattering 
along 6* (right scale), (b) and (c) Double logarithmic repre- 
sentations of the intensity and the width, respectively, versus 
the reduced temperature. 



that stabilizes the charge density modulation in ZrTca 
and measured its temperature dependence. It is an 
acoustic a* polarized phonon branch with a mostly trans- 
verse dispersion along c*. The Kohn anomaly is ob- 
served already at high temperatures of 3-4 Tp and only 
for a sharp momentum space region around qp. On ap- 
proaching the phase transition it becomes giant and the 
phonon freezes into a static modulation of the lattice. 
The temperature dependence of the softening does not 
follow mean field predictions. The analysis of superstruc- 
ture diffraction intensity and diffuse scattering correla- 
tion lengths lead to the conclusion that the phase transi- 
tion is continuous, but governed by an order parameter 
of dimension n = 2. 
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